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Abstract 

It has been suggested that the QPO phenomenon in LMXB's could be explained 
when the central compact object is a strange star. In this work we investigate 
within a standard model for disk accretion whether the observed clustering of spin 
frequencies in a narrow band is in accordance with this hypothesis. We show that 
frequency clustering occurs for accreting strange stars when typical values of the 
parameters of magnetic field initial strength and decay time, accretion rate are 
chosen. In contrast to hybrid star accretion no mass clustering effect is found. 
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1 Introduction 



The recent discovery of new phenomena in the physics of X-ray binaries is 
giving a new and powerful tool to probe dense hadronic matter ( [Miller ct| 
al., 1998[ ) and possibly to infer the existence of a deconfmed phase of quark 
matter in the cores of the central accretors in these systems. Various signals 
of a phase transition to deconfmed quark matter have been suggested in the 
form of peculiar changes of observables, such as the pulse timing ( |Glendenning 
|fc Weber, 2001 ; Poghosyan et al., 2001 ) , and the thermal evolution of isolated 
pulsars QSchaab et al., 1997| ; [Page et al., 2000 E| ; [Blaschke et al., 200 1| ). Even 
more intriguing than the existence of a quark core in a neutron star, is the 
possible existence of a new family of compact stars consisting completely of 
a deconfmed mixture of up (it), down (d), and strange (s) quarks, together 
with an appropriate number of electrons to guarantee electrical neutrality. In 
the literature, such compact stars have been referred to as strange quark stars 
or shortly strange stars (SS), and their constituent matter as strange quark 
matter (SQM). In a series of recent papers ( pey et al., 1998| ; [Li et al., 1999] ) 
it has been argued that the compact objects in some X-ray sources are likely 
strange star candidates. Particularly interesting candidates are the compact 
stars in the newly discovered millisecond X-ray pulsar SAX J1808. 4-3658, and 
in the atoll source 4U 1728-34. 



In the present paper, we model the spin evolution of an accreting strange 
star in a binary stellar system, from the so called "death line" up to the mil- 
lisecond pulsar phase. We explore the dependence of the SS spin evolution 
upon the mass accretion rate and upon those physical quantities which regu- 
late the temporal evolution of the torque acting on the spinning star. We try 
to constrain these quantities to have a population clustering of strange stars 
in agreement with the spin frequency distribution for observed Z sources in 
LMXBs with kHz QPOs. Possible signatures to distinguish strange stars from 
ordinary neutron stars are also briefly discussed. 



2 Accretion model and magnetic field evolution 



We consider the spin evolution of a strange star under mass accretion from 
a low-mass companion star as a sequence of stationary states of configura- 
tions (points) in the phase diagram spanned by the angular velocity Q and 
the baryon number N. The process is governed by the change in angular mo- 
mentum j(n, n) = i(n, n) n 

^-(j(N,n)) = K cxt , (i) 
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where I(N, Q) is the moment of inertia of the star and 



K ext = VGMM 2 r - N out (2) 



is the external torque due to both the specific angular momentum transferred 
by the accreting plasma and the magnetic plus viscous stress given by iV ou t = 
KH 2 r~ 3 , k — 1/3 QLipunov, 1992[ ). For a star with radius R and magnetic field 



strength B, the magnetic moment is given by fi = R 3 B. The co-rotating 
radius r c = (GM/Q 2 ) is very large (r c ^> ro) for slow rotators. The inner 
radius of the accretion disc is 



r 



R , ii < fj,. 
0.52 r A , /j, > /i. 



c 



where \i c is that value of the magnetic moment of the star for which the disc 
would touch the star surface. The characteristic Alfven radius for spherical 

1 /T 

accretion with the rate M = mN is = ^2/i" 4 GMM 2 j ( |Bhattacharya| 



fc van den Heuvel, 1991|) . Since we are interested in the case of fast rotation 
for which the spin-up torque due to the accreting plasma in Eq. (2) is partly 
compensated by N out , eventually leading to a saturation of the spin-up, we ne- 
glect the spin-up torque in N ont which can be important only for slow rotators 
(|Ghosh fc Lamb, 1979| ). 



From Eqs. (1), (2) one can obtain a first order differential equation for the 
evolution of angular velocity 

dQ K cxt (N, Q)-K int (N,Q) 

dt i(N,n) + n(di(N,n)/an) N ' u 

where the internal torque term defined as 

tf MWn ) =n *«°n . (4) 



V dN 



Solutions of (3) are trajectories in the Q — N plane describing the spin evo- 
lution of accreting compact stars . Since for the hybrid stars I(N, Q) exhibits 
characteristic functional dependences ( [Blaschke et al., 200 1| ) at the deconfine- 
ment phase transition line iV cr it(f2) we expect observable consequences in the 
P — P plane when this line is crossed. In the case of SS we have no expecta- 
tion of the mass clustering due to absence of phase transition effects in strange 
matter. 
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In our model calculations we assume that both the mass accretion and the 
angular momentum transfer processes are slow enough to justify the assump- 
tion of quasistationary rigid rotation without convection. For a more detailed 



description of the method and analytic results we refer to |Chubarian et al 
(2000| ) and the works of [Hartle fc Thorne (1968] ), as well as [Sedrakian fc 
Chubarian (19681) . 



The time dependence of the baryon number for the constant accretion rate N 
is given by 



N(t)=N(t ) + (t-t )N . 



(5) 



For the magnetic field of the accretors we consider the exponential decay 
( Bhattacharya fc van den Heuvel, 1991 ) 

B(t) = [B(0) - B^] exp{-t/r B ) + #oo • (6) 



We solve the equation for the spin-up evolution (3) of the accreting star for 
decay times 10 7 < Ts[yr] < 10 9 and initial magnetic fields in the range 0.2 < 
5(0) [TG] < 4.0. The remnant magnetic field is chosen to be B^ = 10" 4 TGQ 
( |Pageet ai, gD0Da| ). 

At high rotation frequency, both the angular momentum transfer from accret- 
ing matter and the influence of magnetic fields can be small enough to let the 
evolution of angular velocity be determined by the dependence of the moment 
of inertia on the baryon number, i.e. on the total mass. This case is similar 
to the one with negligible magnetic field considered in (|Shapiro fc TeukolskyJ 



1983| ; [Burderi et al., 199li| ; |(Jmibarian et al., 200U| ), where /i < fj, c in Eq. (3), 



so that only the so called internal torque term (4) remains. 



3 Equation of state for strange quark matter 



To describe the properties of strange quark matter, we used a recent model for 
the equation of state (EoS) derived by Pey et al. (199"8 ). This model is based on 
a dynamical density-dependent approach to confinement. This EoS has asymp- 
totic freedom built in, shows confinement at zero baryon density, and decon- 
finement at high density. In this model the quark interactions is described by 
a color- Debye-screened interquark potential originating from gluon exchange, 
and by a density-dependent scalar potential which restores chiral symmetry at 
high density (in the limit of massless quarks). This density-dependent scalar 

1 1 TG= 10 12 G 
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Table 1 

Ground state properties of SQM for the equations of state used in the present work. 
(E/A)g S (in MeV) is the energy per baryon, p gs (xlO 14 g/cm 3 ) the mass density, 
and rig S (fm -3 ) the baryon number density. The two equations of state SSI and SS2 
differ for the choice of the parameter v entering in the expression of the in-medium 
quark masses. 

EoS (E/A) gs Pgs ^gs 



SSI 
SS2 
B6O150 



926 
836 



12.3 0.779 
14.1 0.858 
4.6 0.295 



potential arises from the density dependence of the in-medium effective quark 
masses M q , which in the model of Pey et aL (1998 ) are taken to depend on 
the baryon number density tie according to 



M q (n B ) =m q + 310 (MeV) sech(u-), 



(7) 



where q(= u,d,s) is the quark flavor index, no = 0.16 fm -3 is the normal 
nuclear matter density, and v is a parameter. The effective quark mass M q (ns) 
goes from its constituent mass value at zero density, to its current mass m q 
as ub goes to infinity. Here we consider two different parameterizations of the 
EoS, which correspond to a different choice for the parameter v. The equation 
of state SSI (SS2) corresponds to v — 0.333 (y = 0.286). These two models for 
the EoS give absolutely stable SQM according to the strange matter hypothesis 
QBodmer, 197Ij |Witten, 1984] ), see Tab. 1. 



Medium dependent mechanisms for confinement and their consequences for the 
EoS of quark matter, have been explored by many authors using different QCD 
motivated phenomenological models ( [Blaschke et al., 1990| , |1999| ; piaschke &] 
dandy, 2000 ; Drago et al., 1996' ). In addition to the previous model for the EoS, 
we make use of the MIT bag model EoS ( [Farhi fc Jaffe, 1984]) for strange quark 
matter for non-interacting quarks, with strange quark mass m s = 150 MeV, 
massless u and d quarks, and with bag constant B = 60 MeV/fm 3 (hereafter 
the B6O150 EoS). 



4 Results and discussion 



In Fig. 1 we show evolutionary paths of accreting strange stars in the mass- 
radius (MR) plane (top panels) corresponding phase diagrams (lower panels) 
for two different EsoS (SSI right panels and B60i 50 left panels). In each panel 
we show two trajectories of a strange star initially rotating with frequency 
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Q(0) = 0.001 Hz; for which the initial magnetic field is B(0) = 7 TG and 
its decay time is tb = 10 s yr. The mass accretion rate onto the star is M = 
lO~ 9 M /yr. The solid lines show the evolution of strange star configurations 
with initial gravitational mass M(0) = 1.4M , the dashed lines show that 
of configurations with initial baryon mass iV(0) = 1AN Q . In the insets of 
Fig. 1 we show the stable branches for configurations rotating with maximal 
frequency (dash-dotted lines) in compare to static ones (dotted lines). 




Fig. 1. Evolutionary paths for strange stars in the mass-radius plane (top panels) 
and in the frequency-baryon number plane (lower panels) for the equations of state 
SSI (left panels) and B6O150 (right panels), see text. 

From the evolution paths Q(N) one can see that in all cases of initial masses 
the magnetic braking force is strong enough to stop fast spin-up of the star 
and saturate the frequency of rotation. This can lead to an effect of frequency 
clustering. However such effect for strange stars is correlated with the accretion 
rate. 
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In Fig. 2, we show the result of our calculation for the spin evolution of the 
accreting strange star with EoS SSI in the P — P diagram, where P = 2n/Q 
is the period of rotation. The parameters of the accretion model are chosen 
such as to correspond to values extracted from observations made on LMXBs, 
which are divided into Z sources with M ~ lCr 8 M Q /yr and A(toll) sources 
with M ~ 10~ 10 Af Q /yr (|Bhattacharya fc van den Heuvel, 1991| ; |van der Klis 
Mlendenning fe Weber, 200l| ). 
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Fig. 2. P - P diagrams for accreting strange stars with equation of state SSI. We 
show the dependence on the initial magnetic field -B(O) (upper panel), mass accretion 
rate M (lower panel), and magnetic field decay time tb (middle panel). 

In Fig. 2 we explore the sensitivity of the model to changes of the parameters 
with respect to the set used in Fig. 1 (B(0) = 7 TG, tb = 10 8 yr, M = 
lO~ 8 M /yr) for accretors with initial baryon mass N(0) = 1.4 TVq. In the 
upper panel we vary the inital magnetic field B(0) = 5, 7, 8 TG, in the middle 
one the magnetic field decay time tb = 10 7 , 10 s , 10 9 yr and in the lower panel 
the accretion rate M = 10~ u , 1(T 9 , l(T 8 M Q /yr. 

We see from this figure that changes of the accretion parameters not only 
shift (choice of the initial magnetic field) and deform (choice of the magnetic 
field decay time) the interval where the spin- up is saturated (dip in P), but 
for some cases this effect can be washed out by a variation of the accretion 
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rate. This phenomenon means that the existence of a frequency clustering for 
strange stars requires a strong limitation of possible values of the accretion 
parameters. 

In Fig. 3 we plot the 'Waiting time' r = \P/P\ ( [Poghosyan et al., 200 1| ) of the 
strange star as a function of the spin frequency, for the EsoS SSI and B6O150. In 
order to obtain an enhanced waiting time (frequency clustering) in the interval 
220 Hz < v < 380 Hz which corresponds to recent observations (|Glendenning| 
fc Weber, 200 1| ), we have to choose the following parameters: B(0) = 2.5TG, 



r B = 10 8 yr, M = lO- 9 M /yr for SSI model and 5(0) = 3 TG, r B = 10 s 
yr, M = lO _8 M /yr for B6O150 model. For both cases the initial gravitational 
mass is the same 1.4M , which corresponds to 1.83iV initial baryon mass 
for SSI and 1.71iV for B6O150 model. The main difference between both of 
these scenarios for the frequency clustering is the initial baryon mass and the 
mass accretion rate. An independent determination of these quantities to a 
sufficient accuracy could thus rule out one of the compact star models. At 
present, this could be done only for a minority of objects ( [Lamb fc Miller 



2001] ) and bears some model dependence. ( |Poghosyan et al., 200 1| ). 
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Fig. 3. Waiting times as a function of spin frequency for spin-up evolution of strange 
stars with equations of state B6O150 an d SSI. 
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5 Conclusions 



We have investigated the question whether the clustering of spin frequencies 
which has been observed for the compact objects in LMXBs is consistent with 
the hypothesis that at least some of these objects are strange stars. We have 
applied a standard model of magnetic disc accretion and find that population 
clustering in a narrow band of spin frequencies between 220 < z/[Hz] < 380 
can occur for typical parameter values of the model. Inspection of the response 
to parameter variations shows that the lowering of the magnetic field decay 
time can wash out the effect as well as a change in the accretion rate. The 
changes in the initial magnetic field leave the waiting time distribution rather 
unchanged but shift the interval of period clustering. 

On the other hand, the observation of frequency clustering alone is no indi- 
cation for the presence of a strange star. A standard bag model EoS shows a 
similar waiting time pattern as the strange quark matter EoS does. For the 
EsoS investigated in this paper no clustering of masses has been obtained. 
This is a striking difference to hybrid star configurations ( [Poghosyan et aE 



200 1| ), which could be used in order to constrain further our approaches to 
the EoS of superdense matter provided a sufficiently large sample of LMXBs 
could be observed and their frequencies and masses could both be extracted 
with sufficient accuracy. 
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